function y=nansem(x)
y=nanstd(x) ./ sqrt(sum(~isnan(x)));
end